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ABSTRACT 



As a contribution to quantum optics in the vicinity of surfaces we study the single 
atom spontaneous emission in a linear chain of two-level atoms. The electromagnetic field 
is thereby treated with the help of integro-differential equations which take into account 
the interaction with the other atoms in the chain. The life time of the excited atom, the 
frequency shift of the atomic transition and the angular distribution of emitted photons are 
worked out. They depend on the position of the emitting atom. As compared with the single 
atom in free space, considerable modifications occur for atoms a few interatomic distances 
away from the ends of the chain. 
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1. INTRODUCTION 



During the last decades there has been a growing interest in the behaviour of atoms 
and their spontaneous emission in particular near the surfaces of media. See, for example, 
the reviews [TTI] — jTJ. The common method to solve boundary problems in quantum optics 
is based on the solution of Maxwell equations for the macroscopic electric field strength. 
Thereby the respective boundary conditions must be specified. The spatial geometry of the 
medium and the optical properties are given by the dielectric permittivity e. The initially 
excited atom is characterized by the transition frequency and the transition dipole moment. 

At present physical situations are in the center of the interest in which optical phenomena 
take place in the nearest vicinity of the media interfaces. At distances much less than the 
wave length of the spontaneous photons, various surface effects play an important role. For 
example, the influence of evanescent modes on an atom near a vacuum-dielectric interface 
has been studied in PH. The atomic spectral characteristics and the spatial distribution of 
photons emitted spontaneously by the atom located in vacuum at small distances from the 
interface differ in principle from corresponding characteristics when the atom is located at 
large distances from the surface. In the theoretical analysis carried out in |8],|9| the standard 
approach was used: the electromagnetic field was described by the macroscopic Maxwell 
equations and the medium was considered to be a continuous system of radiators. 

The excistence of evanescent modes is not the only surface effect which may modify 
influence on the process of spontaneous emission. A different type of influence shows up if 
the discrete atomic structure of the medium is taken into account. In this case it can be 
demonstrated that even in a spatially uniform dielectric medium the optical properties of 
the surface layers differ significantly from the bulk properties. This is due to the violation of 
the translational symmetry at the surface. For a classical and a quantum optical treatment 
10]- Jl2| and Jl^JT^I , respectively. Here the characteristics of the spontaneous emission 



sec 



of an atom for the case of semi-infinite medium were obtained. 
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Regarding the description of the medium, in the papers ||13| , p^| a very promising mi- 
croscopic approach to the electromagnetic field was introduced which is based on integro- 
differential equations for the local field strength which is generated by a system of single 
radiators (atoms, molecules) situated at specific places with particular physical proper- 
ties [|T^,|nj. In the discussion below we will continue to apply the local field approach to the 
photon field but change to a medium with different geometry. 

In recent years the investigation of the optical properties of linear molecular aggregates 
has obtained a lot of attention. In many cases such systems may be successfully modeled as 
linear chains of two-level atoms interacting by dipole-dipole interactions. At present the non- 
linear optical properties of these systems ]17l[18] and the process of collective spontaneous 



emission (superradiation) [|19[ are sufficiently well understood. The results in papers |T8|JT9 



are again obtained on the basis of the classical electromagnetic theory. They indicate that 
in the radiation field the evolution of the atoms, located very close to the ends of the chain, 
is essentially different from the corresponding evolution of atoms located far away from the 
ends. This motivated the investigations of the present paper. 

Below we study in the framework of the Weisskopf-Wigner approximation and based on 
the local field approach the process of spontaneous emission of a two-level atom embedded 
in a linear chain of iV — 1 ground state atoms. We will show that the spectral characteristics 
of the atom, as the life time of the excited state and the frequency shift of the atomic 
transition (Lamb shift), depend strongly on the atomic position if the atom is located very 
close to the ends of the chain while this is not the case at a larger distance. The angular 
distribution of the spontaneously emitted photons is also investigated. We thereby rely 
on the mathematical formalism developed in papers [|13|,|20|| but in contradistinction to the 



works [18||n| the electromagnetic field will be quantized. 

The paper is organized as follows: In Section 2 the general formalism is presented which 
is based on the Heisenberg equations for the atomic operators and the field operators and 
an integro-differential equation for the local electromagnetic field strength. In Section 3 
we calculate the electromagnetic mode functions for the linear atomic chain. Section 4 is 



devoted to the investigation of the life time and energy shift of an excited atom embedded 
in the linear chain. In Section 5 we investigate the angular distribution of the spontaneously 
emitted photons. Concluding remarks are made in Section 6. 



2. THE FUNDAMENTAL EQUATIONS 

We firstly describe a system of N two-level atoms located at arbitrary positions r/ 
(I = 1, N). We assume the atom with number j (J = 1, N) with transition frequency 
ujq to be in the excited state while the other atoms with transition frequency ui ~ ujq are in 
their ground states. Our aim is to study the process of spontaneous emission emission of 
this atom. As compared with the free space case (one excited atom in vacuum) this process 
will be modified by the presence of the other atoms which are assumed to be polar izable and 
may therefore have an induced dipole moment. The initially excited atom which performs 
a dipole transition is during this process in a dipole-dipole interaction with all the other 
atoms. In our approach this will be represented by the adequately chosen vector potential 
operator. We then have for the initially excited atom at the position Tj the following effective 
Hamiltonian [BI|,B2,0: 



H = -huj a 3 + — rf u d A(r j , t)a 2 + ^ ^fcC^ A c kA , (2.1) 

1 kA 

where 02, cx 3 are the Pauli matrices, = do/ do with do being the matrix element of 
the electric dipole moment operator is calculated with the aid of the unperturbed wave 
functions. c^ A , CkA are the photon creation and annihilation operators for the kA-mode. We 
assume that the interaction of the atom with the photon field isn't too strong, so that the 
Weisskopf-Wigner approximation holds. 

The vector potential operator A(rj,t) needs a particular attention. At the position of 
the j-th atom it may be written in the Heisenberg representation in the form of the mode 
decomposition 



2ixhc 2 

A(r J ,t)=Y^g' kX i kX (r J )c kX (t) + h.c., g' kX = ] J^-—, (2.2) 



where Vr is a quantization volume and fkA is the mode function for the physical situation 
sketched above. It corresponds to the wave vector k and the polarization index A with 

ujk = ch- 
in order to obtain the mode functions it is a common approach to solve Maxwell equations 
with appropriate boundary conditions. This procedure has the disadvantage that apart 
from evanescent waves || various surface effects are not taken into account. Because we are 
interested in the following in such effects which go back to the violation of the translational 
symmetry at the boundary of the medium, we choose a different approach which takes into 
account the physical properties of single atoms located at specific places. This method 
turns out to be particularly suited for the study of quantum optical effects which happen 
very close to the boundary. Instead of the macroscopic Maxwell equations we use thereby 
an integro-differential equation for the local classical field strength [|15|,[n| . In the case of the 



discrete system of polarizable atoms it has the form of a system of coupled linear equations 

fiaOj) = e A exp(ifcsrj) +^ rot rat i a I f kA (rj)G ! (.Ry) , j = l,...,N, (2.3) 

with G(R) = exp(ikR)/R and Rij = \ri — Tj\. Differentiation in (2.3) is carried out with 
respect to coordinates Tj. ati is the direction independent polarizability of the /-th atom. e A 
is the unit vector of polarization, s = h/k is the unit vector in the direction k of the mode 
in question and k is the wave number. Note that because of the second term in (2.3) which 
represents the influence of the other atoms (I ^ j) at the position r,- of the excited atom, the 
modulus | fkA (r?) I will depend on rj. This will lead to the boundary effects which show up 
in modified characteristics of the process of spontaneous emission as compared with the free 
space case (only the first term in (2.3)). The solution of equations (2.3) will be discussed 
below. 

Regarding the physical interpretation of (2.3) it has been shown in that the quantum 
analogue of this integro-differential field equation may be obtained within the framework 
of third-order QED effects. Thereby the evolution of an iV-atom system in the radiation 
field has been investigated in the dipole approximation. Each atom can emit or absorb one 



real photon independently of the others. In third order perturbation theory this process 
includes the interaction between the atoms through the exchange of virtual photons. The 
latter results in the second term of (2.3). 

To work out the life time of the excited atom and the frequency shift of its atomic 
transition, we turn to the Heisenberg equations. Introducing the notations 



Aj = Au d = J29kx(rj)c k x + h.c. , gyxfe) = 9kX f kx(rj)n d 

kA 



(2.4) 



we obtain as in 23 



uJodo 

(7i = -uqo-2 + 2——Ajcr 3 



cr 2 = lo (?i , 



(2.5) 



0"3 



he 



CkA 



.uj d c 



he 



Equations (2.5) agree structurally with those for the spontaneous emission of a single 
atom in the free space. The solution for a not too short time t after the beginning of the 
spontaneous emission is 



Ck\(t) = c{ r x ee {t) + i-f—gt\{Yj) [(y + {t)C{uk + uj ) - a(t)C{uJk - w ) 
(a 3 (t)) = -1 + [1 + (a 3 (0))] exp [-t/rfo)] , 
(a + (t)) = (a+(0)> exp [i (u + Sfa)) t) exp H/2r(r,)] , 



(2.6a) 



with 



1 



Afr, 



2tt 



/ u d 
{ he 

( ^od[ 



El^kA 
/ kA 

/ kA 



(r^)! 2 [S(u k - u ) - 5(u k + u )] 



P 



(2.66) 



he ) ^ '"—v j" \ UJk -\- U j Uk — ojQ, 
Here cr + , a are the raising and lowering operators of the two- level atom, respectively; ((x) 



is a Heitler function [23|. Brackets (■ ■ •) denote vacuum averaging over the state \^/) = 
\vacuum) |$), where |<&) is an arbitrary state of two-level atom and \vacuum) is the state 
without photons. 
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Regarding the interpretation of (2.6a) we see that the photon operator refers in its first 
term to the free space radiation field and in the additional term to the photon source part. 
The latter field is distorted by the medium (atoms with I ^ j) in accordance with the 
coordinate dependence of the mode functions fk\(rj). Furthermore, it can directly be read 
off from (2.6a) that r is the life time of the excited state and A the frequency shift of the 
atomic transition. Both now depend on the position of the excited atom so that surface 
effects will become evident. 



3. MODE FUNCTIONS OF A ONE-DIMENSIONAL CHAIN OF ATOMS 



We specialize the situation described above and consider a one-dimensional straight chain 
of N atoms with equal polarizability a located with interatomic distance ao on the negative 
part of the z-axis, so that r/ = (0, 0, zt) with z\ = — a (/ — 1) and I — 1, • • • , N. Performing the 
differentiations in (2.3) we obtain the following system of equations for the mode functions: 

fj = e A exp [-iks z a (j - 1)] + (3.1) 

v f. 2 fj-|n, ,, 3/fn 2 -f, 3ffn,-f, ) .... 
f£[ a \l-j\ a 2 (l-j) 2 ag|/-j| 3 J 

where = fkA( r /) = (ft \ fi \ ft)- n z ^ s a un ^ vec tor in the direction of the 2;-axis. It is 
convenient to introduce spherical coordinates. The direction of the wave vector k of the 
mode becomes: 

s = (sin 6* cosy?, sin^ sin ip, cos 9) , (3.2) 

where 9 is an incident angle which is counted from the positive direction of the axis z 
(0 < 9 < 7r), tp is azimuth angle (0 < 9 < 2ir). Then polarization vectors will have the form: 

= (sin(/9, — cosy?, 0) , ey = (— cos 9 cos <p, — cos^ sin ip, sin 9) . (3.3) 

Because of the symmetry of the problem the result of calculations will not depend on p and 
we may restrict to p = n: 

s = (-sin#,0,cos#) , ej_ = (0,1,0) , e M = (cos^, 0, sin 9) . (3.4) 
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Projecting in (3.1) on the coordinate axes, we then obtain the following system of coupled 
scalar equations: 

f] I- j ■ ^(7,17 , (3.5a) 

1+3 

1+3 

with 

"0 



exp (ikao\l — j\) , (3.56) 

where Ej = exp [—ika (j — 1) cos#] corresponds to free field. 

The equations (3.5) can be interpreted as follows: The atom j undergoes the spontaneous 
emission. The influence on this by an other atom I ^ j in the chain is given by the factors 
Cij of (3.5b). The problem contains three lengthscales: the interatomic distance a , the 
wavelength A = 2ix/k of the mode in question and the polarizability a 1 / 3 . The physics is 
characterized by two dimensionless parameters a/a^ and a k. The parameter ka is assumed 
to be small because we are dealing with optical quantum transitions. The common factor 
a/al gives the strength of the interaction between the atoms. The phase of the different 
C[j is determined by the wavelength of the mode according to the distance a \l — j\ of the 
contributing atom /. The term proportional to \l — j\~ 3 represent the Coulomb part of the 
interaction. It falls of quickly with the distance of the atom I so that retardation effects do 
not play a role. The other terms are the retardation terms. With ka® 1 their contribution 
will be small. The overall picture is therefore that the sums in (3.5a) will essentially reduce 
to only a few terms. We want to discuss this in detail. 



k Qjq 



ika 



exp (ika \l — j\) 



C? 3 = 



2a 



ikao 



+ 



3.1. Numerical results 
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To be able to judge the approximations we are going to made below, we make a numerical 



approach first. Equations (3.5) are solved with the iterative method of Zeidel P5] . For 
the wavelength A = 628 nm (corresponding to the resonance frequency Uq) we choose the 
interatomic distance ao to be ao = A/628 so that ka$ <C 1 is fulfilled. To make the resulting 
figures easily interpretable we choose as number of atoms in the chain iV = 629 so that its 
length becomes L = A. 

The resulting y- component of the mode functions is shown in Fig.l for different values 
of the normalized polarizability C = a/a^. They demonstrate that there is no visible 
modification of the wave length of electromagnetic radiation. This means that the chain 
behaves like a dielectric medium with a refractive index equal to unity. The reason for the 
fact that phase properties of the field don't undergo perceptible changes is the following: 
According to the exponential factor in (3.5b) essential phase shifts at the position of the j-th 
atom in comparison with the phase of the free field may occur only due to the interactions of 
the j-th atom with far distant atoms, i.e., if the quantity \l—j\ is sufficiently large. However, 
for large \l — j\ the coefficients C\j tend to zero because of the square brackets in (3.5b). 
This fact explains the discovered property of the chain. 

Because the amplitudes of the real and imaginary parts of the field in Fig. 1 a-c are 
equal to one another, the direction of propagation of the field in a one-dimensional chain 
coincides with the direction of propagation of the external field. This fact may apparently 
be considered as a consequence of the property discussed above. 

We add a remark: The refractive index forms as a result of a summation of phase shifts 
produced by every sufficiently close atom. Therefore, owing to the character of the dipole- 
dipole interaction, there may well be a non-trivial refractive index for systems where the 
atoms are not arranged as a linear chain as discussed here. For higher dimensional arrays 
every atom has a larger number of neighbours. A more detailed discussion of this problem 
will be given in a subsequent paper. 

The physical situation described above changes though a few interatomic distances away 
from the ends of the chain. There we observe according to Fig. Id rapid oscillations going 



back to the violation of the translational symmetry. The amplitudes of these oscillations and 
the distances from the ends of the chain at which they occur depend on the polarizability a. 
With the increase of a the distances, at which boundary effects are essential, increase too. 

To represent these oscillations which happen on top of the behaviour in the free space 
case (a = 0, first term in (3.5a)) we make the ansatz 

f? = T?E jCO Be , f] = T]E 3 , (3.6) 

f* = T*E jS me , 7/ 7/ , j = l,...,N. 

The amplitude factors T? and Tj 2 remain to be specified, j denotes the position of emit- 
ting atom. Analytical expressions can be obtained on the basis of physically simplifying 
assumptions. We turn to this now. 

3.2. Nearest neighbours approximation 



If a is small enough we have shown that it is justified to take into account in the coupled 
equations (3.5) only the interactions with one or two nearest neighbours. This leads to 
analytical approximate solutions of (3.5) which will later be compared with the numerical 
results. 

Restricting first to the interaction with only one nearest neighbour, equations (3.5) take 
the form: 

f? = E jCOS 9 + C? (//_! + , 

f? = E j + CP 1 (fV_ 1 + fy +1 ) , (3.7) 
/; = ^sin0 + C*(i-_ 1 + /* K ) , 

W — °i - --3 > u i - -j - -^W • 
Here we have taken into account only the Coulomb interaction of the dipoles because of 
ka <C 1 in the optical range. We seek for the solution of the system (3.7) in the form 

f» = T x Ej cos6 , fy = T%Ej , f; = T*E 3 sin6 , (3.8a) 

10 



for the positions 2 < j < N — 1 and 

f? = T?E jC os9 , f? = T?E j , f^T^Ej sin 9 , (3.86) 

for the ends of the chain (j = l,N). Substituting (3.8a), (3.8b) into (3.7) and using again 
ka <d 1, it is easy to obtain the unknown amplitudes. The result is 

T S = T »=YT2C ' T ° = ~ 2T ° ' ^ 
1 + c 

rpx rpy r F z ctrpx 

h ~ h ~TT2d ' h ~~ Ul ■ 

To improve the approximation we now take into account interactions with two nearest 
neighbours. In this case the system of equations (3.5) reduces to 



i+2 J ' 



// = E j cos o + ci (/;• , + /;.,) + c| ( jj 2 + /; 

fj = E j + c l (fj-i + + c l (fj-2 + //+ 2 ) , (3.10) 
// = ^ sin + C? + + C| (#_ 2 + f> +2 ) , 



la _ 1 a 

with coefficients Cf, Cf, defined in (3.7). We look for a solution of the system (3.10) in 
the form (3.8a) for 3 < j < N — 2, in the form (3.8b) for j = 1, N, and in the form 

f? = T?E j <x*e , f] = T y 2 E 3 , f* = T*E jS m9 , (3.11) 

at the positions j — 2, N — 1. Substituting (3.8), (3.11) into (3.10) we then obtain 

rpx rpy 1 

°~1 + 2.25C ' 
1 + 0.125C-2.125C 2 

2 2 — 1 2 — 



[1 + 2.25C) (1 - C 2 ) ' 



T x_ T y_ l 1-125 + 0.125C-2C 2 

A " A - 1 " 6 (1 + 2.25C)(1-C 2 ) ' {6A2) 
7? = -27? , Z = 0,1,2 . 



11 



The result for atoms at the ends of the chain and far from the ends (T^ y ' z ) is given 
in Table 1 for different values of the normalized polarizability C = a/a^. Note that the 
disagreement between the numerical calculation and the calculation based on the nearest 
neighbour approximation is only of the order of a few percent. This approximation reflects 
therefore very well the underlying physical processes. 



3.3. The infinite dipole chain 



To complete the discussion we treat the infinite dipole chain. The limit N — > oo corre- 
sponds to the neglect of boundary effects. To obtain an analytical solutions of the system 
(3.5) we make the "ansatz" (3.8a) for every j. Substituting (3.8a) into (3.5) and omit- 
ting small terms of the order ka and higher, we obtain the following expressions for the 
amplitudes 

1 2 

T o = T o = i + 2CC(3) ' T ° 2 = "l + 2CC(3) ' (3 - 13) 

where £(3) = 1.202 is the third order Riemann ^-function. The results T of (3.13) agree 
exactly with the numerical results of the last coulomn of Table 1 which refer to the inner 
atoms. Thus, due to predominant character of the Coulomb part of the dipole-dipole inter- 
action, boundary effects in a linear chain don't influence the behaviour of the mode functions 
far from the ends of the chain. 



4. LIFE TIME OF THE EXCITED STATE AND FREQUENCY SHIFT OF THE 

ATOMIC TRANSITION 

We are now prepared to work out the life time of the excited state and the frequency 
shift of the atomic transition (Lamb shift) in a linear chain of atoms in making use of formula 
(2.6) and the results obtained in section 3. Here we shall not limit the treatment to nearest- 
neighbour approximation but use the general expression (3.6). Taking into account (3.3), 
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the mode functions may be written as 



(4.1) 



fkAfo) = E\j exp (iksTj) , 

E\\j = (-7J cos 6 cos ip, -Tf cos sin ip, T* sin 0) , 
E±j = 7J(sin</?, -cos</?,0) , j = l,---,N . 
We substitute (4.1) into (2.6) and replace the mode summation by a summation over polar- 



izations and an integration over wave vectors according to 

V R 

A 



£ 

kA 



(2ttc) 3 

where is a solid angle element in the direction of emission of spontaneous photon. Inte- 
grating over the frequencies we get 



T 



I\) Ti 8lT V J4tt 



(4.2) 



where ri is the life time of a single atom in free space. The unit vector u d of (2.1) which 
defines the orientation of the transition dipole moment of the emitting atom at rj may be 
written as 

u d = (sin d cos ip d , sin d sin ip d , cos d ) . (4.3) 

Substituting (4.1), (4.3) into (4.2) and integrating over the angles we finally obtain for the 
position dependent life time r{vj) 



{Tf)\m 2 6 d +(T*)\os 2 6 d 



(4.4) 



t(tj) n 

Corresponding calculations for the frequency shift of the atomic transition give an analogous 
result 

A(r,) = A 1 [(lj) 2 sin 2 d + (t*) 2 cos 2 d ] , (4.5) 

where Ai is the free space frequency shift. This is the intended result. For the different T 
see (3.12) or Table 1. 

The dependence of the life time r{vj) on the number j of the initially excited atom is 
shown in Fig. 2 for the case d = n/2. We see that the life time of the excited state may 
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be considerably larger for atoms near the end of the chain as compared with an atom in 
free space. Only for atoms several interatomic distances away from the end, the value of the 
life time in an infinite chain is reached which is still larger than the one in free space. All 
these effects increase with increasing C = a/a^. Note that r -1 (rj) and A(i-j) depend in the 
same characteristic way on the angle 9d which the transition dipole moment forms with the 
direction of the linear chain. 



5. ANGULAR DISTRIBUTION OF THE SPONTANEOUSLY EMITTED 

PHOTONS 

Finally we investigate the angular distribution of the photons spontaneously emitted 
by an atom in a linear chain of atoms. To do so we introduce the operator of the photon 
number in the kA-mode: n^x = c kA c kA- Going back to section 2 and the Hamiltonian (2.1) 
and calculating the corresponding commutators, we obtain the Heisenberg equation which 
governs the time evolution of the operator of the number of photons emitted by the j-th 
atom 

n k \(r j: t) = iZ^- [gk\(rj)ckx(t) - ^ A ( r j)4 A (t)] . (5.1) 

The photon creation and annihilation operators are defined by the formula (2.6). After 
elementary transformations the vacuum mean value of the time derivative of the photon 
number operator may be wtitten in the form 

he 



(n kX (r J ,t)) = 2nh^\ |, kA (r,)| 2 x 



x 



(cr + a) j S (uj k - luq) + (<7<7 + ) j 5 (u k + u ) . (5.2) 



We study as usual the intensity in the direction of emission s (3.2) at the time t for the 
atom located at the position rf 

V -I' 

H r 3,s,t) = j^^J2j^k{nkx(r j ,t))ujlduj k . (5.3) 
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Substituting (5.2) into (5.3) we get an expression which is similar to the one obtained in f23 



I(Tj,S,t) 



0"0 i + 



2nc 3 



(5.4) 



The fk A( r j) are still unspecified. Therefore, the formula (5.4) is the general expression for 
the emission intensity and it may be used for the investigation of the angular distribution 
of spontaneously emitted photons in media with arbitrary atomic distribution. In the case 
of a single atom in free space we have fj^r^) = e^exp (iksjrj) and the formula (5.4) gives 
the well-known result for the intensity of dipole radiation: 



(5.5) 



To discuss the linear chain, we define a new unit vector 

-V2 , 



cos 2 9 d 



rpz 

3 



sin 9 d cos ip d , sin 9 d sin ip d , ^ cos 9 d 



(5.6) 



Note that the direction of depends on IJ/TL? and therefore on the position of the atom 
in the chain. Making use of this notation and of the formulas (3.3) and (4.1), we obtain in 
this case finally 



As a result it follows from (5.7) that the characteristic structure of the angular distri- 
bution of the photons emitted by the atom embedded in a linear chain is analogous to the 
one of a single atom in free space [equation (5.5)]. The dipole-dipole interaction leads to 
two modifications: There is a factor in (5.7) which diminishes the intensity uniformly 

for all directions. Secondly the well-known free-space angular distribution is rotated by an 
angle 7 given by 

cos7 = u d u^ . (5.8) 

In the optical range we have T? = — 2T? so that the angle 7 becomes the same for all atoms. 
Substituting (4.3), (5.6) into (5.8) we get 

1 - 3 cos 2 9 d 



2 r 



1 



(5.7) 



cos 7 



VI + 3 cos 2 9 d 
15 



(5.9) 



The dependence of the rotation angle 7 on the orientation of the transition dipole moment 
is shown in Fig. 3. 

6. CONCLUDING REMARKS 

We have studied the single-atom spontaneous emission in a linear chain of two-level 
atoms. Particular attention has been payed to emitting atoms situated very close to the 
ends of the chain. It is shown that the respective life time, Lamb shift and the angular 
distribution of the intensity of the spontaneously emitted photons change considerably as 
compared with atoms in the middle of the chain or atoms in free space. 

To treat the quantized electromagnetic field we did not rely on the macroscopic Maxwell 
equations with boundary conditions and space dependent refractive index, but used the 
integro-differential equations for the local field instead. From the point of view of QED the 
respective integro-differential equations for the local field take into account that the emission 
of a real photon may be accompagnied with an exchange of virtual photons with the other 
atoms of the chain. The discrete structure of the medium and the physical properties of 
the single atoms at given positions are thereby naturally represented. The fundamental 
interatomic interaction is the dipole-dipole interaction. We remark that a chain of dipoles 
also appears in the problem of quantum optics in the three-dimensional microscopic cavities, 
compare [p6jJT^| . However, the behaviour of these image dipoles does not agree with the one 
of the real dipoles which build up the medium in our case. 
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FIGURE CAPTIONS 



Fig.l. Numerical results for the position dependence of the mode function in a linear 
atomic chain extending from z = to z = —628 nm. The interatomic distance ao is 1 nm. 
C = a/al with the atomic polarizability a. The wavelength of the free space photon field is 
assumed to be equal to the length of the chain. It can be read off from the Fig.a-c that in a 
chain the resulting wavelength is the same. There are strong oscillations a few interatomic 
distances away from the ends of the chain which are shown in detail in Fig. Id. The solid 
and the dashed line correspond to the real and imaginary part, respectively. 

Fig. 2. Life time r of a single excited atom in a linear chain, j denotes the number of 
the position of the atom counting from the end of the chain. T\ is the life time in free space. 
For the atomic transition the wave length A = 628 nm and the orientation of transition 
dipole moment was assumed to be perpendicular to the direction of the chain (0 d = ir/2). 
The other parameters are chosen as for Fig.l. 

Fig. 3. For an atom in a chain the free space distribution of emitted photons is rotated 
by an angle 7 which depends on the orientation of the transition dipole moment (angle 6a) 
relative to the direction of the chain. The corresponding dependence is shown in this figure. 



19 



Table 1. Amplitude factors of the mode functions obtained by approximations and 
merical calculations. 





interaction with one 
nearest neighbor 


interaction with two 
nearest neighbors 


numerical solution 
of the system (3.5) 


c 


Tf 


rpx 
1 2 


rpx 

1 


Tf 


rpx 

1 2 


rpx 

J o 


Tf 


rpx 

1 2 


rpx 

J o 


0.1 


0.917 




0.833 


0.908 


0.817 


0.816 


0.902 


0.812 


0.806 


0.2 


0.857 




0.714 


0.846 


0.675 


0.690 


0.839 


0.668 


0.675 


0.3 


0.813 




0.625 


0.807 


0.555 


0.597 


0.800 


0.546 


0.581 


0.4 


0.778 




0.556 


0.786 


0.445 


0.526 


0.787 


0.425 


0.510 


0.5 


0.750 




0.500 


0.784 


0.333 


0.471 


0.817 


0.262 


0.454 
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